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Abstract We report recent progress in the problem of distinguishing convex subsets 
of cyclotomic model sets A by (discrete parallel) X-rays in prescribed A-directions. 
It turns out that for any of these model sets A there exists a 'magic number' ota such 
that any two convex subsets of A can be distinguished by their X-rays in any set of 
ota prescribed A-directions. In particular, for pentagonal, octagonal, decagonal and 
dodecagonal model sets, the least possible numbers are in that very order 11,9, 11 
and 13. 



1 Introduction 

Discrete tomography is concerned with the inverse problem of retrieving informa- 
tion about some finite set in Euclidean space from (generally noisy) information 
about its slices. One important problem is the unique reconstruction of a finite point 
set in Euclidean 3-space from its (discrete parallel) X-rays in a small number of 
directions, where the X-ray of the finite set in a certain direction is the line sum 
function giving the number of points in the set on each line parallel to this direction. 

The interest in the discrete tomography of planar Delone sets A with long-range 
order is motivated by the requirement in materials science for the unique recon- 
struction of solid state materials like quasicrystals slice by slice from their images 
under quantitative high resolution transmission electron microscopy (HRTEM). In 
fact, in |[T2l . lfT4l a technique is described, which can, for certain crystals, effec- 
tively measure the number of atoms lying on densely occupied columns. Clearly, 
the aforementioned density condition forces us to consider only A-directions, i.e. 
directions parallel to lines through two different points of A. 
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In the quasicrystallographic setting, the positions to be determined form a finite 
subset of a nonperiodic Delone set with long-range order (more precisely, a model 
set Q, iTPTl ). Model sets possess a dimensional hierarchy, i.e. they allow a slicing 
into planar model sets. In fact, many of the model sets that describe real quasicrys- 
tallographic structures allow a slicing such that each slice is an n-cyclotomic model 
set, the latter being (planar) Delone sets contained in the additive subgroup of the 
Euclidean plane generated by the nth roots of unity; cf. J6], 0, and lfl5l for 
details. It therefore suffices to study the discrete tomography of these cyclotomic 
model sets. In practice, the cases n = 5, 8, 12 are of particular interest. In the present 
text, we shall mainly focus on the larger class of cyclotomic Delone sets. 

Since different finite subsets of an n-cyclotomic model set A may have the same 
X-rays in several directions, one is naturally interested in conditions to be imposed 
on the set of directions together with restrictions on the possible finite subsets of A 
such that the latter phenomenon cannot occur. Here, we consider the convex subsets 
of A and summarise recent results in the problem of distinguishing convex subsets 
of A by X-rays in prescribed A-directions. It turns that there are four prescribed A- 
directions such that any two convex subsets of A can be distinguished by the corre- 
sponding X-rays, whereas less than four A-directions never suffice for this purpose. 
Much more novel is the result obtained in collaboration with M. SpieB that there is 
a finite number m^ such that any two convex subsets of A can be distinguished by 
their X-rays in any set of m^ prescribed A-directions. Moreover, the least possible 
numbers m^ in the case of the practically most relevant examples of n-cyclotomic 
model sets A with n = 5, 8 and 12 only depend on n and are in that very order 11,9 
and 13. This extends a well-known result of R. J. Gardner and P. Gritzmann J5) on 
the corresponding problem for the crystallographic cases n = 3,4 of the triangular 
resp. square lattice A (with least possible number ota = 7 in both cases) to cases that 
are relevant in quasicrystallography. 

The intention of this text is to provide an easy to read guide to the results of iflOl 
with a view towards practical applications. Detailed proofs, related results and an 
extensive list of references can be found there. For the algorithmic reconstruction 
problem in the quasicrystallographic setting, the reader is referred to 0, Q. 



2 Cyclotomic Delone sets 

Throughout the text, the Euclidean plane is identified with the complex numbers. For 
Z € C, z denotes the complex conjugate of z. Further, we denote by K\ the smallest 
subfield of C that contains the rational numbers as well as the union of A — A and its 
image A — A under complex conjugation. Recall that A is called a Delone set if it is 
both uniformly discrete and relatively dense. For n G N, we always let = e 2m ' n , 
a primitive nth root of unity in C. Then, the smallest subfield of C that contains the 
rational numbers as well as is the nth cyclotomic field denoted by Q(C„). The 
latter ist just the Q-span of the nth roots of unity and thus contains the Z-span Z[£„] 
of the nth roots of unity. Recall that a homothety of the complex plane is given by 
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z i-4 Az + 1 , where X <E R is positive and t S C. For the purpose of this text, the 
following rather abstract definition provides a convenient framework. 

Definition 1. Let n > 3. A Delone set A C C is called an n-cyclotomic Delone set if 
it satisfies the following properties: 

(n-Cyc) K\ is contained in Q(Cn) ■ 
(Horn) For any finite subset F of K\, there is a homothety h of 
the complex plane that maps the elements of F to A . 

Further, A is called a cyclotomic Delone set if it is an n-cyclotomic Delone set for a 
suitable n > 3. 

Standard examples of n-cyclotomic Delone sets are the n-cyclotomic model sets, 
which were also called cyclotomic model sets with underlying 1,-module Z[£„] in |9l 
Section 4.5] and are defined via the canonical cut and project scheme that is given 
by the Minkowski representation of the Z-module Z[£„]; see ll9l [lOl for details. 
These sets are certain Delone subsets of the Z-module Z[£„] and range from periodic 
examples like the fourfold square lattice (n ~ 4) or the sixfold triangular lattice 
(n = 3) to nonperiodic examples like the vertex set of the tenfold Tubingen triangle 
tiling (n = 5), the eightfold Ammann-Beenker tiling (n = 8) or the twelvefold shield 
tiling (n = 12); see 03 Figure 1], |9j Figure 2] and FigureQ]below for illustrations. 
Note that the vertex sets of the famous Penrose tilings of the plane fail to be 5- 
cyclotomic model sets but can still be seen to be 5 -cyclotomic Delone sets; see |fl] 
and references therein. 



3 Determination of convex subsets by X-rays 

Let (t\,t2,t?,,t4) be an ordered tuple of four distinct elements of RU {°°}. Then, its 
cross ratio (t \ , t 2 , h , t^) is the nonzero real number defined by 

{h ^ hM) ~ {t,-t 2 )(u- h y 

with the usual conventions if one of the f, equals °°. 

The unit circle in C is denoted by S 1 and its elements are also called directions. 
For a nonzero complex number z, we denote by sl(z) the slope of z, i.e. sl(z) = 
—i(z -z)/(z + z) e R U {°°}. Let A be a subset of C. A direction u e S 1 is called 
a A-direction if it is parallel to a nonzero element of the difference set A — A = 
{v — w | v, w G A} of A. By construction, the cross ratio of slopes of four pairwise 
nonparallel A-directions is an element of the field K\ (~l R. In case of n-cyclotomic 
Delone sets A, these cross ratios are thus elements of the field Q(£„) nR. 

Definition 2. • Let F be a finite subset of C, let u € S be a direction, and let 
be the set of lines in the complex plane in direction u. Then the ( discrete parallel) 
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X-ray of F in direction u is the function X U F : — >• No = N U {0}, defined by 



• Let & be a collection of finite subsets of C and let 17 C § be a finite set of 
directions. We say that the elements of & are determined by the X-rays in the 
directions of U if different elements of & cannot have the same X-rays in the 
directions of U, 

One can easily see that no finite set of pairwise nonparallel A-directions suffices 
in order to determine the whole class of finite subsets of A by the corresponding 
X-rays (9J. It is therefore necessary to impose some restriction on the finite subsets 
of A to be determined. It has proven most fruitful to focus on the convex subsets 
of cyclotomic Delone sets. The latter are bounded (and thus finite) subsets C of A 
satisfying the equation C = conv(C) n A, where conv(C) denotes the convex hull of 
C. One has the following fundamental result which shows that one has to choose 
the set U of A-directions in such a way that certain convex polygons cannot exist; 
cf. |9j Proposition 4.6 and Lemma 4.5]. Here, for a finite set U C S 1 of directions, 
a nondegenerate convex polygon PcCis called a U -polygon if it has the property 
that whenever v is a vertex of P and u€U, the line in the complex plane in direction 
u which passes through v also meets another vertex v' of P. P is called a V '-polygon 
in A, if its vertices lie in A. Note that the proof of direction (ii)=^(i) needs property 
(Horn) and see FigureQ]for an illustration of the other direction (i)=>(ii). 

Theorem 1. Let A be a cyclotomic Delone set and let U C S 1 be a set of two or 
more pairwise nonparallel A-directions. The following statements are equivalent: 

( i) The convex subsets of A are determined by the X-rays in the directions of U. 

( ii) There is no U -polygon in A 

In addition, i/card(f/) < 4, then there is a U -polygon in A. 

Employing Darboux's theorem on second midpoint polygons J4) together with 
a blend of sophisticated methods from the theory of cyclotomic fields and previous 
results obtained by Gardner and Gritzmann 0, one obtains the following deep result 
on U -polygons; cf. iflOl Theorem 5.7]. Note that the proof heavily relies on property 
(n-Cyc). 

Theorem 2. Let n>3 and let A be an n-cyclotomic Delone set. Further, let (/cS 1 
be a set of four or more pairwise nonparallel A-directions and suppose the existence 
of a U -polygon. Then the cross ratio of slopes of any four directions ofU, arranged 
in order of increasing angle with the positive real axis, is an element of the finite set 
^lcm(2n,l2)(Q(Cn) Hi) of numbers in the field Q(£„) flR that can be written in the 



X U F(£) ^card(Fn^). 



form 
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Fig. 1 The boundary of a U -polygon in the vertex set A of the twelvefold shield tiling, where U is 
the set of twelve pairwise nonparallel A-directions given by the edges and diagonals of the central 
regular dodecagon. The vertices of A in the interior of the [/-polygon together with the vertices 
indicated by the black and grey dots, respectively, give two different convex subsets of A with the 
same X-rays in the directions of U. 



where (^1,^2,^3,^4) is an element of the set 

{(kifoMM) € N 4 |fc 3 < < 1(2 < A4 < lcm(2«, 12) — 1 and k\ +1(2 = £3 +£4} • 

Moreover, card(t/) is bounded above by a finite number b„ E N that only depends 
on n. In particular, one can choose £>3 = b\ = 6, b$ = 10, b% = 8 andb\2 = 12. 

Theorems Q] and |2]now immediately imply our main result on the determination 
of convex subsets of cyclotomic Delone sets; cf. iflOi Theorem 5.11]. 

Theorem 3. Let n > 3 and let A be an n-cyclotomic Delone set. 

(a) There are sets of four pairwise nonparallel A-directions such that the convex 
subsets of A are determined by the corresponding X-rays. In addition, less than 
four pairwise nonparallel A-directions never suffice for this purpose. 

(b) There is a finite number m„ £ N that only depends on n such that the convex 
subsets of A are determined by the X-rays in any set of m n pairwise nonparallel 
A-directions. In particular, one can choose m^ = m<± = 1, = 11, m% = 9 and 
OT12 = 13. 



By Theorems Q] and |2]above, it suffices for Part (a) to take any set of four pair- 
wise nonparallel A-directions such that the cross ratio of their slopes, arranged 
in order of increasing angle with the positive real axis, is not an element of the 
finite set <r ^icm(2«,i2)( ( Q(C;i) HR); cf. [10, Corollary 4.10] for concrete results 
in the practically most important cases n = 5,8, 12 of quasicrystallography. 
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Our above analysis heavily relies on the assumption of ideal data and is therefore 
only a very first step towards a satisfactory tool for materials science. Further, it 
would certainly be interesting to abandon the slice by slice approach and work, 
for a Delone set A in Euclidean 3-space, with A-directions in general position in- 
stead; compare the approach to 3D reconstruction of atomic arrangements presented 
in ifTTl . In that case, it might well be that seven is a universal magic number for the 
determination of convex subsets by X-rays; cf. @ Problem 2.1]. 
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